Majorana modes in time-reversal invariant s-wave topological superconductors 
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We present a time-reversal invariant s-wave superconductor supporting Majorana edge modes. 
The multi-band character ol the model together with spin-orbit coupling are key to realizing such a 
topological superconductor. We characterize the topological phase diagram by using a partial Chern 
number sum, and show that the latter is physically related to the parity of the fermion number of the 
time-reversal invariant modes. By taking the self-consistency constraint on the s-wave pairing gap 
into account, we also establish the possibility of a direct topological superconductor-to-topological 
insulator quantum phase transition. 
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Since Majorana suggested the possibility for a fermion 
to coincide with its own antiparticle back in 1937 [1], the 
search for the Majorana particle has catalized intense ef- 
fort across particle and condensed-matter physicists alike 
Q. Particles either constitute the building blocks of a 
fundamental physical theory or may effectively emerge 
as the result of the interactions of a theory. A Majorana 
fermion is no exception to this principle, with neutrinos 
potentially epitomizing the first view [3j, and localized 
quasiparticle excitations in matter illustrating the second 
[4| . Remarkably, Majorana fermions can give rise to the 
emergence of non-Abelian braiding Q. Thus, in addition 
to their significance for fundamental quantum physics, 
interest in realizing and controlling Majorana fermions 
has been fueled in recent years by the prospect of imple- 
menting fault-tolerant topological quantum computation 
@j 13 • As a result, a race is underway to conclusively 
detect and characterize these elusive particles. 

A variety of condensed-matter systems hosting local- 
ized Majorana elementary excitations have been pro- 
posed, notably certain quantum Hall states @ and so- 
called topological superconductors (TSs) [1, 01 ■ Unfortu- 
nately, these exotic states of matter require the explicit 
breaking of time-reversal (TR) symmetry and their phys- 
ical realization seems to be at odds with existent mate- 
rials. Such is the case, for instance, of superconductors 
with p x + ip y spin-triplet pairing symmetry. This has not 
prevented researchers to pursue creative proposals that 
rely on a combination of carefully crafted materials and 
devices. Fu and Kane in particular, suggested the 
use of a (topologically trivial) s-wave superconducting 
film on top of a three-dimensional topological insulator 
(TI), which by proximity effect transforms the non-trivial 
surface state of the TI into a localized Majorana excita- 
tion 0, 11 1 (see also 12 1 for related early contributions). 
While experimental realization of this idea awaits further 
progress in material science, alternative routes are be- 
ing actively sought, including schemes based on metallic 
thin-film microstructures, quantum nanowires, and semi- 
conductor quantum wells coupled to either a ferromag- 



netic insulator, or to a magnetic field in materials with 
strong spin-orbit (SO) coupling (l3| . 

Our motivation in this work is to explore whether a 
path to TSs exists based on conventional bulk s-wave 
spin-singlet pairing superconductivity. We answer this 
question by explicitly constructing a model which, to the 
best of our knowledge, provides the first example of a 
2D TS with s-wave pairing symmetry, and supports Ma- 
jorana edge modes without breaking TR symmetry [lij ]. 
The key physical insight is the multi-band character of the 
model, in the same spirit of two-gap superconductors [H| , 
but with the SO coupling playing a crucial role in turn- 
ing a topologically trivial two-gap superconductor into a 
non- trivial one. Our results advance existing approaches 
in several ways. First, multi-band systems clearly expand 
the catalog of TI and TS materials. Following the dis- 
covery of s-wave two-band superconductivity in MgB2 
in 2001, a number of two-gap superconductors ranging 
from high-temperature cuprates to heavy-fermion and 
iron-based superconductors have already been character- 
ized in the laboratory [lij ]. giving hope for a near- future 
material implementation. Furthermore, from a theoreti- 
cal standpoint, our TR-invariant model also supports a 
direct TI-to-TS (first-order) quantum phase transition 
(QPT), allowing one to probe these novel topological 
phases and their surface states by suitably tuning con- 
trol parameters in the same physical system. 

Exact solution with periodic boundary conditions. — 
We consider a TR-invariant two-band Hamiltonian of the 
form H = f/cd + H so + H SVJ + H.c, where 
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j,u=x,y 



(i) 



represent the two-band (c and d) dynamics, the SO in- 
teraction, and s-wave superconducting fluctuations, re- 
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spectively. In the above equations, fi is the chemical 
potential, u c d represents an onsite spin-independent "hy- 
bridization term" between the two bands, fermionic cre- 
ation operators at lattice site j (unit vectors x, y) and 
spin a =t, 4- are specified as c t a or eft , depending on 
the band, and (A c , Ad) denote the mean-field s-wave 
pairing gaps. By letting = (c jt f,Cj i i,dj ! f,d jt i) T , the 
Pauli matrices t v and a v act on the orbital and spin part, 
respectively. Notice that we have implicitly assumed that 
the intraband SO coupling strengths obey A c = — A^ = A. 
In this way, in the limit where [i = = A c = Ad, H re- 
duces to a known model for a TI [l7j |. 

For general parameter values and periodic bound- 
ary conditions (PBC), H can be block-diagonalized 
by Fourier transformation in both x and y. That 
is, we can rewrite H = h X)k(^iL^kAc — 4/i), with 

= ( C k,t ' C k,i' ^k,f> ^KV C -k,t' C -k,^' ^-k,f ^-k,J.)' anC ^ 

i?k an 8 x 8 matrix. An analytical solution exists in 
the limit where the pairing gaps are 7r-shifted, A c = 
= A, since decouples into two 4x4 matri- 
By introducing new canonical fermion operators, 

= 7f ( C k,o- + dk.cr), kk,<T = ^(Ck.o- - <4:,<r), We 
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-k,f> "-k,4.J u -k^' ^-k,f a k,^ ^k,t 

H 2 k , with iJj k , i?2 k being TR of one another 



and F k = 



iA(Jy 

-m^a z + n + Ai 



Here, Ak = — 2A(sinfc z: ,sinfcy), mt = — 2t(cosk x + 
cosfc y ), and c? = (a x ,a v ). The excitation spectrum ob- 
tained from diagonalizing either H[ k or H' 2 k is 



e«,k — ± 



-fi 2 + |A k | 2 ±2 



m k £7 2 - 



^ 2 |Ak| 2 , (2) 



where the order ei.k < £2,k < < 63. k < £4,k is as- 
sumed and Q, 2 = n 2 + A 2 . QPTs occur when the gap 
closes (e 2j k = 0, for general A^O), leading to the crit- 
ical lines determined by — ±f2, where the criti- 
cal modes k c € {(0, 0), (0, 7r), (7r, 0), (n, tt)}. It is worth 
noticing that through a suitable unitary transformation 
(see Eq. (4) of Ref. [Hi), the SO interaction in Eq. flTJ) 
is formally mapped into p x + ip y and p x — ip y intraband 
interaction, hinting at the existence of non-trivial topo- 
logical phases, as we demonstrate next. 

Topological response. — Since H preserves TR invari- 
ance, bands which form TR-pairs have opposite bulk 
Chern numbers (CNs) C n , leading to £ neoccupied c n = 
(including both H [ k and H' 2 k ). Thus, introducing a new 
Z2 topological invariant is necessary in order to distin- 
guish between trivial and TS phases. In Ref. [lij], the 
parity of the sum of the positive CNs was considered, 
whereas in Ref. [2(]j an integral of the Berry curvature 
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FIG. 1: (Color online) Topological characterization of the 
phase structure of Hamiltonian H via the partial CN sum C+ 
as a function of u c d and A, with t = 1 and arbitrary A / 0, 
for representative chemical potentials [i — (top) and /i = — 1 
(bottom). The black (dashed) line represents an insulator or 
metal phase, depending on the filling, with A = 0. CNs are 
calculated for (N x ,N y ) = (100, 100) lattice sites. Note that 
we may have two pairs of edge modes with C+ = 0. 



over half the Brillouin zone for all the occupied bands was 
used. Here, we propose a different Z2 invariant which is 
guaranteed to work in the presence of TR: taking advan- 
tage of the decoupled structure between TR-pairs, we use 
the CNs of the two occupied negative bands of H[ k only 
(say, C\ and C2) and define the following parity invariant: 



P c = (-i) m ° d ^+) ; C+ = d + C 2 



(3) 



Let |?/vk) denote the band-n eigenvector of H[ k . Then 
the required CNs C n , n — 1,2, can be computed as [21] 



" 



dk y lm(d kj: ipn,k\dk y ipn,k)- ( 4 ) 



The resulting topological phase structure is shown in 
Fig. [1] by treating the pairing gap A as a free control 
parameter. In an actual physical system, A cannot be 
changed at will, but only be found self-consistently by 
minimizing the free energy (or ground-state energy at 
zero temperature). While we shall return on this issue 
later, we first focus on understanding the physical mean- 
ing of the above invariant and on establishing a bulk- 
boundary correspondence for our model. 

Interestingly, there is a direct connection between 
the invariant Pc and the fermion number parity of the 
TR-invariant modes. Without loss of generality, let 
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fj, = 0, and focus on the ground-state fermion number 
parity of the four TR-invariant points in the first 
Brillouin zone, k c . Since H[ k and H' 2 k are decoupled, 
we need only concentrate on the ground-state parity 
property of H[ k . Let us introduce the new basis: 

i a k c tl VaC )' ^k c J VaC )' l vaC )' a k c t^L il vaC )l- ^ n tmS 

sis, H[ k becomes -ffi,k r = mk c a z ®Aa x , with eigenvalues 
±mk c , ±A, and an identical matrix for H' 2 k in the TR- 

,t' ° 



basis 
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vac) 



~ a k c ,^k c ,t 



| vac}}. 
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When |m,k e | > |A|, the ground state of each 
mode k c is in the sector with odd fermion parity, 



P kc = e "H Q kc ,T a k c,T+ h kc ,4-°kc,4.> = _x, otherwise it is 
in the sector with even fermion parity, Pk c = 1- 
By analyzing the relation between |?7ikJ and |A| for 
each k c , we can see that the TS (trivial) phases with 
Pc — —1(1) correspond to the ground state with 
rik c ^ k c = Pf = —1(1)- Thus, our Z2 invariant 
coincides with the fermion number parity of the four 
TR-invariant modes from one representative of each 
Kramer's pairs, consistent with the fact that only a 
partial CN sum can detect TS phases in the presence of 
TR symmetry. While the relation between non-trivial 
topological signatures (such as the fractional Josephson 
effect) and the local fermion parity of Majorana edge 
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states has been discussed in the literature [g, [22 
invoking the fermion number parity of the TR-invariant 
modes in bulk periodic systems to characterize TS phases 
has not, to the best of our knowledge. 

Open boundary conditions and edge states. — A hall- 
mark of a TS is the presence of an odd number of pairs of 
gapless helical edge states, satisfying Majorana fermion 
statistics. Thus, in order to understand the relation be- 
tween Pc (or Pp) and the parity of the number of edge 
states, i.e., a bulk-boundary correspondence, we study 
the Hamiltonian H on a cylinder. That is, we retain 
PBC only along x, and correspondingly obtain the exci- 
tation spectrum, e„.fe x , by applying a Fourier transforma- 
tion in the a;-direction only. For simplicity, let us again 
focus on the case /1 = 0. The resulting excitation spec- 
trum is depicted in Fig. [2] for representative parameter 
choices. Specifically, for odd Pc (C+ = 1 in panel (a) 
and C+ = — 1 in panel (b), respectively), H supports one 
TR-pair of helical edge states on each boundary, corre- 
sponding to the Dirac points k x — (a) and k x = ir (b). 
Different possibilities arise for even Pc- While C+ = 
can clearly also indicate the absence of edge states, in 
panel (c) one TR-pair of helical edge states exists on 
each boundary for both Dirac points k x — 0,7r (for a to- 
tal of two pairs, as also explicitly indicated in Fig. [I). In 
panel (d) (C+ = 2), both TR-pairs of helical edge states 
correspond to the Dirac point k x = instead. Since, as 
remarked, our Hamiltonian exhibits particle-hole symme- 
try, the equation j e = J- Cn k holds for each eigen- 
value Cn.kxi where j e is the associated quasi-particle 
annihilation operator. Thus, for zero-energy edge states, 




FIG. 2: (Color online) Excitation spectrum of Hamiltonian 
H on a cylinder for fj, — 0, t = 1, A = 1. Panel (a), C+ = 1: 
A = 2, u cd = 3; Panel (b), C+ = -1: A = 2.5, u cd = 2; Panel 
(c), C+ = 0: A = 2, u cd = 1; Panel (d), C+ = 2: A = 0.8, 
u c d = 1.5. Note that the bulk gap scales as min(A,A). The 
number of lattice sites (N x ,N y ) = (40, 100). 



7 = 7q, indicating that the edge states in our system 
satisfy Majorana fermion statistics. 

Phase diagram with self- consistent pairing gap. — 
Within BCS mean-field theory, let V = V^.k' > de- 
note the effective attraction strength in each band. Then 
the pairing gap A = A c = -F(c k ,tC-k4> = -A d , 
and the ground-state energy can be written as E g = 
2N x N y (A 2 /V) + Ek( £ i,k + e 2 ,k - 2//). The first (con- 
stant) term is the condensation energy, which was ne- 
glected in H. By using Eq. @ and minimizing E g , we 
obtain the stable self-consistent pairing gap A as a func- 
tion of the remaining control parameters [24j . The result- 
ing zero-temperature phase diagram is shown in Fig. [3] 
For fi — (top panel), the average fermion number is 
consistent with half-filling, and thus with an insulating 
phase when A = 0. In particular, when < \u c( i\ < 4, 
the ground state is known to correspond to a TI phase 
171 ]. Interestingly, without self-consistency, the TI can- 
not be turned into a TS directly, as shown in the top panel 
of Fig. [T] However, after self-consistency is taken into 
account, the topologically trivial phase with C+ = ±2 
disappears, and a first-order QPT can connect the two 
phases. For fj, = — 1 (bottom panel), the average fermion 
number is found to be less than half-filling, realizing a 
metallic phase when A = 0. Derivatives of the ground- 
state energy indicate that all QPTs, except the TI-to-TS 
phase transition, are continuous. 

Discussion. — A number of remarks are in order. First, 
while the choice of SO coupling strengths and s-wave 
pairing gaps obeying A c = — A<j and A c = — A^ affords a 
fully analytical treatment, relaxing these conditions may 
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FIG. 3: (Color online) Phase diagram as a function of u c d 
and V with the pairing gap A calculated self-consistently. 
The magnitude of A is represented by a color, whose scale is 
indicated on the side. The number of lattice sites (N x , N y ) = 
(80,80). 

be necessary to make contact with real materials. Nu- 
merical results on a cylinder show that the level crossing 
of the Majorana edge states in the TS phase is robust 
against perturbations around A c = — A^, including the 
possibility that the SO coupling vanishes in one of the 
bands. TS behavior also persists if |A C | — \Ad\ ^ 0, as 
long as the phase difference between pairing gaps is n. 
In the presence of a phase mismatch e, edge modes are 
found to become gapped, with a minimal gap that scales 
linearly with e. Interestingly, however, preliminary re- 
sults indicate that adding a suitable Zeeman field can 
allow (at the expense of breaking TR invariance) gapless 
Majorana excitations to be restored, with a precise tun- 
ing of the phase difference being no longer required. It is 
also worth noting that one can reinterpret the band index 
in H as a layer index, and so H may be thought of as de- 
scribing a bilayer of superconductors with phase-shifted 
pairing gaps, and an interlayer coupling H c a. Beside es- 
tablishing a formal similarity with the scenario discussed 
by Fu and Kane [loj . such an interpretation may offer 
additional implementation flexibility, as the possibility 
to control the superconducting and SO couplings by an 
applied gate voltage has been demonstrated recently (2|| ■ 
Second, we have thus far restricted to 2D systems in or- 



der to simplify calculations. Preliminary results indicate 
that a qualitatively similar behavior (that is, the possi- 
bility of even/odd numbers of pairs of gapless Majorana 
surface states) also exists for 3D systems obtained from 
a natural extension of our 2D Hamiltonian. It is espe- 
cially suggestive to note that a it phase shift in the order 
parameter across two bands is also believed to play a key 
role in iron pnictide superconductors [26{ , hinting at pos- 
sible relationships between TS behavior and so-called s± 
pairing symmetry. While a more detailed investigation is 
underway, it is our hope that multi-band superconduc- 
tivity may point to new experimentally viable venues for 
exploring topological phases and their exotic excitations. 
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